Correlated vs Ferromagnetic State in Repulsively Interacting Two-Component Fermi 

Gases 
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Whether a spin- 1/2 Fermi gas will become ferromagnetic as the strength of repulsive interaction 
increases is a long-standing controversial issue. Recently this problem has been studied experimen- 
tally by Jo et al, Science, 325, 1521 (2009) in which the authors claim a ferromagnetic transition is 
observed. This work is to point out the results of this experiment can not distinguish whether the 
system is in a ferromagnetic state or in a non-magnetic but strongly short-range correlated state. A 
conclusive experimental demonstration of ferromagnetism relies on the observation of ferromagnetic 
domains. 



Itinerant ferromagnetism is a common phenomenon in 
nature, but not yet well understood. Rigorous examples 
of itinerant ferromagnetic ground state have only been 
obtained for a few specific cases. For instance, Nagaoka 
shows that for infinite strong repulsive interaction, in a 
bipartite lattice the ground state is ferromagnetic if one 
hole is doped into a half-filled system [l[ . Lieb shows for 
a half filled bipartite lattice, the ground state of repul- 
sively interacting fermions has non-zero spin if the num- 
ber of total lattice site of each sub-lattice is not equal 
Q. Mielke [3| and Tasaki [3] propose a class of mod- 
els whose single particle ground states have degeneracy, 
and show they become ferromagnetic with repulsive in- 
teractions. However, there is no conclusive results for a 
generic dispersion and filling number. 

Stoncr considered spin-1/2 fermions with short range 
interactions, spin polarization can lower the interaction 
energy since two spin align fermions will not interact due 
to the Pauli exclusion principle, while it costs the ki- 
netic energy. With Hatree-Fock approximation, one can 
conclude that when UAf(E-p) > 1 there exists a second- 
order ferromagnetic phase transition where U is the 
interaction strength and J\f(E F ) is the density-of-state 
nearby the Fermi surface. This is known as Stoner crite- 
ria. For s-wave scattering, this condition corresponds to 
kpa s > 7r/2, where k-p is the Fermi momentum, and o s is 
the s-wave scattering length. Higher order perturbation 
of interactions will lower the critical value of fcpa s , and 
may change the transition to first order 

Many authors have proposed to study itinerant ferro- 
magnetism transition using two-component Fermi gases 
where a s can be tuned by Feshbach resonance [|| . Based 
on the physical picture above, in a trapped system one 
should observe non-monotonic dependence of the kinetic 
energy with the increase of a s , namely, the kinetic energy 
shall first decrease before ferromagnetic transition due to 
the expansion of the cloud, and then increase after the 
transition. The inelastic collision rate shall first increase 
and then decrease as different components begin to sep- 
arate spatially Recently, a beautiful experiment by 
Jo et al have observed all these monotonic features, 
and the agreement between experiment and ferromag- 



netic theory [ll| leads to the claim that this has shown 
experimentally a ferromagnetic transition in continuum 
without particular requirement of lattice and band struc- 
ture 0. 

However, the itinerant ferromagnetic issue is in fact 
more complicated than this. The question is, whether 
spin polarization is the only way to reduce interaction en- 
ergy. The answer is no. In the content of Hubbard model, 
Gutzwiller constructed his famous projected wave func- 
tion as Yiii^ ~ 7 7 n iT n «l)l^ / o}) where \^o) is free fermion 
Fermi sea, and i is the index of the lattice site. The pro- 
jection operator Yii(^ ~ V n i]' n il) (v > 0) suppresses the 
probability of having two fermions at the same lattice 
site, and consequently reduces on-site interaction energy 
This state is non-magnetic if j^o) is chosen as non- 
magnetic state. Hereafter we shall call this state "corre- 
lated state" to distinguish it from "ferromagnetic state" . 
Nevertheless, we shall note this state is not an exotic 
state but still a Fermi liquid state, we use the term "cor- 
related state" in the sense that the projection operator 
introduces strong short-range correlation into this state. 
In continuum, a Jastrow factor can play the role of the 
projection operator. 

In the Hubbard model, using the projected wave func- 
tion as a variational wave function, Gutzwiller shows that 
at low-density, the correlated state has lower energy than 
a ferromagnetic state 12J. An alternative view is that 



the short-range correlation, which has been ignored in 
the Hatree-Fock and perturbation treatment, will signif- 
icantly renormalize down the interaction. Kanamori ar- 
gued that the up-bound of the effective interaction should 
correspond to the kinetic cost to put a node in the wave- 
function where two fermions overlap, which should alway 
be finite even when bare interaction goes to infinite, and 
he also argued that the renormalized interaction is not 
sufficient for ferromagnetic transition at low density [l3| , 
which is supported by some later calculations 

In short, the key of the itinerant ferromagnetism prob- 
lem is whether the system will choose spin polarization 
or building up short-range correlation to reduce interac- 
tion energy as the strength of interaction increases. The 
advantage of cold atom is to provide an opportunity for 
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a direct quantum simulation of the Stoner model, and 
hopefully can settle the issue of itinerant ferromagnetism 
experimentally. So the question comes to whether the ex- 
periment of Ref. [HI has conclusively settled the issue. 
The answer is no. The purpose of this Rapid Communi- 
cation is to point out a non-magnetic "correlated" state 
can explain the main observation of Ref. [Io| equally 
well as a ferromagnetic state, in another word, from the 
existing experimental results, it is very hard to distin- 
guish whether the system is in a "correlated" state or in 
a ferromagnetic state. Further experimental efforts are 
required to distinguish them. 

Equation-of-state for a "correlated" state. Let us first 
consider two-component fermions in free space (without 
optical lattice and harmonic trap), the Hamiltonian is 
given by 



H 



E 



2m 



(1) 



where v(rj — rj - ) is a short-range pairwise interacting 
potential. For a non-polarized free Fermi sea j^o) = 

Det(e lkir i r )Det(e lkjr j ), the kinetic energy of each compo- 
nent is given by E® in = 3£?pn./5, where Ep = /i 2 fc|/(2m), 
n is the density of each component, and n = kp/(6ir 2 ). 
For a Fermi sea, the interaction energy is proportional to 
the k = Fourier component of v(r) (denoted by vq), i.e. 



£? nt = (*o| $>(rj -rj)|*„) 



v n 2 . 



(2) 



Away from a Feshbach resonance, E° nt = Airh 2 a s n 2 /m. 

Now we consider Gutzwiller's projected wave function 
in continuum = V\^o) as a class of varational states. 
With the projection operator, the probability of having 
two spin-opposite fermions closely changes from n 2 to 
(1 — g)n 2 , and the interaction energy decreases if g > 
and increases if g < 0, thus the interaction energy shall 
linearly depend on the "projection strength" g as 



4irh 2 a. 
-Bint = n 



(3) 



By dimension analysis the kinetic energy shall be of the 
form 



Skin = <* | 



rf,<7 



-^|*) = \E F n w{9 ), (4) 



where w{g) is a dimensionless function of g. There are 
some simple properties of w{g) one can make use of. For 
g = 0, there is no projection and the free Fermi sea is the 
state that minimizes the kinetic energy, thus w(0) = 2. 
If g 0, both positive and negative g will lead to the 
increase of the kinetic energy, thus g = is the minimum 



of w(g), namely, dw(g)/dg\ g= o = 0. Hence, up to the 
second order of g, one has the form 

E kin = \E F n{2 + ag 2 ) = ^-a n 5/3 (2 + ag 2 
5 m 



(5) 



where a > 0, and a = 3(67r 2 ) 2 / 3 /(407r) = q.36. We shall 
now stress that the purpose of this work is neither to 
rigorously derive this equation-of-state and calculate the 
number of a, nor to prove theoretically that this state 
can energetically do better than a ferromagnetic state. 
Instead, we shall take Eq. [3] and Eq. [S] together as a sim- 
ple "phenomenological " equation-of-state for this class 
correlated state, and the key of work is to point out the 
general behavior of this correlated state in trap, which 
does not depend on the specific value of a, and hereafter 
we shall use a as an unspecified parameter. 

For a given density n and a s , one shall first minimize 
the free energy with respect to g. For ag?! 1 / 3 < 2aao, 
g = a s n 1 / 3 /(2aao). In this regime, 
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FIG. 1: (Color online) The total energy (blue solid line), the 
kinetic energy (red dashed line) and the interaction energy 
(black dotted line) of the "correlated state" as a function of 
a s fcF. For this plot we set a = 1. E^ in is the kinetic energy of 
a free Fermi sea 
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and Ei nt = 0, the total energy 
4ttS 2 



Etnt = 



a n 5/3 {2 + a) 



and the chemical potential 



4irh 2 
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The kinetic, interaction and total energy (in unit of E® in ) 
as a function of a s fcp are illustrated in Fig. [TJ When 
g = 1 at very large fcpa s , the energy of a correlated state 
is lower than a fully polarized ferromagnetic state if a < 
2 5 / 3 -2 ~ 1.17. 

Trapped System. From the discussion above, we have 
obtained the relation /z(n,a s ). For a given a s , one can 
invert this relation to obtain n(/j,,a s ). Considering the 
harmonic trapping potential Vtrap(r) = (muj_(x 2 + y 2 ) + 
muj z z 2 ) /2, we shall use local density approximation to re- 
place /j, with hq — Vtrap(r) and by solving the total number 
of particle constraint / rf 3 rn(/io — V(r), a s ) — N, one can 
obtain /j,q(N, a s ). Then the local fermion density is given 
by n(r) = n(fJ,o(N, a s ) — V(r),a s ). Using the expres- 
sions for kinetic and interaction energy density discussed 




FIG. 2: (Color online) From (a)-(e) is the kinetic energy 
per particle E k = £km/(2iV), the interaction energy per 
particle E\ = £i nt /(2N), the potential energy per particle 
E p = £ P ot/(2AT), the total energy per particle E t = £ to t/{2N) 
and the chemical potential /io as a function of a s ko- fco is the 
Fermi momentum for free Fermi gas at the center of the trap, 
and the energy unit is taken as = fi 2 fco/(2m). (f) the 
three-body loss rate as a function of a s ko- a = 1 for blue 
solid line, a = 0.75 for red dashed line and a = 0.5 for black 
dotted line. 



above, one can compute the total kinetic and interac- 
tion energy as £ kin /int = J d 3 rE kin/int (n(r),a s ), and the 
potential energy is given by £ pot = / d 3 rVt rap (r)n(r), 
and the total energy is £ tot = £kia + £mt + £pot- The 
loss rate is computed in a very phenomenological way as 
r = 2T (a s k a ) e J d 3 rn 3 (r)(l - g(n(r)) [ljj. 

As shown in Fig. [T] for a uniform system the kinetic 
energy for a correlated state monotonically increases for 
any a s > 0. To show whether for small a s the kinetic en- 
ergy will first decrease with the increase of a s in a trapped 
system, we shall note 



kin 



da s 



d 3 r 



<9-Ekin dn d/j, , dE kia 



dn d\i da s 



da s 



(13) 



The first term is negative and the second is positive. It is 
important to note that when a s — > the first term does 
not vanish while the second term does, since /io linearly 
depends on a s while -Ekin quadratically depends on a s , 
therefore the first term is always dominative in small a s , 
which gives d£ k i n /da s < 0, and leads to a non-monotonic 
behavior of kinetic energy. 

We consider the experimental condition as Ref. [Io| . 
i.e. N = 6.5 x 10 5 and u z /uj± = 7/30. The results 
are shown in Fig. [2] Comparing them with the predic- 
tion of a ferromagnetic state, for instance, Fig 1 and 2 
of Ref. [U and Fig. 2 of Ref. dl], they display simi- 
lar non-monotonic behavior and also qualitatively agree 
with the observation of Ref. • This leads to the main 
point of this work, that is, the non-monotonic behavior 
observed in Ref. [Uj is not sufficient to distinguish a fer- 
romagnetic state from a non-magnetic correlated state, 
and thus not conclusive for making the claim of ferromag- 
netic transition. We emphasize that despite of the similar 
non-monotonic behavior, there is no phase transition in 
this scenario. In fact, the suppression of interaction en- 
ergy and the inelastic collision rate due to correlation is 
not surprising in strongly interacting systems. Quantum 
Hall effect and the Tonk gas of one-dimensional bosons 
are two of the examples. Suppression of the three-body 
recombination rate has been observed in one-dimensional 
Bose gas as it approaches the Tonk gas regime [l6| . 

Discussions. There are a few points we would like to 
comment on before ending. First, there are some quanti- 
tative differences between the results of Fig. [5] and that 
from a ferromagnetic theory (for instance, Fig 1 of Ref. 
1). In Fig. [2J the extreme of kinetic energy, potential 
energy and the loss rate are not very close, while they are 
very close in the ferromagnetic theory prediction. And 
there is no maximum in the chemical potential (i.e. cloud 
size) plot of Fig[^e). However, both calculation above 
and the theoretical work of Ref. [1, @, [Tl[ are not quan- 
titatively correct. The important effect of Feshbach res- 
onance and unitary limit of the repulsive interaction is 
not taken into account. For instance, the Hatree-Fock 
energy of a free-Fermi gas is taken as linearly increasing 
with a s , while the accurate Hatree-Fock energy should 
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be smaller and saturates at large a s . The resonance 
physics has to be taken into account seriously for mak- 
ing a quantitative comparison between theory and ex- 
periments, for instance, the value kpa s of kinetic energy 
turning point, and for constructing a correct microscopic 
Fermi liquid theory. And for the correlated state, the cor- 
rection should be treated more seriously rather than the 
phenomenological way presented above, for instance, by 
quantum Monte Carlo simulation. It remains to be seen 
whether these quantitative difference between the pre- 
diction of two scenarios can be used to distinguish these 
two states, when a more careful analysis in the theory is 
done. We leave this for follow up works. 

Secondly, a conclusive experimental evidence of ferro- 
magnetism is the observation of ferromagnetic domains. 
Ref. [Icj fails to observe the ferromagnetic domains. 
They attribute this reason to short lifetime that prevents 
the system to reach equilibrium. However, one should no- 
tice that this system is the same as what has been used to 
study BEC-BCS crossover before. Maybe there is some 
particular physics reason to believe the relaxation time 
is particularly long in this case than in the case of BEC- 
BCS time. If it is the case, the dynamics remains to be 
explored. 
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